
Calhoun 

iniQiuiic^iul Ar{hiv« of tilt Mil vdl Poii^roduiit School 


Calhoun: The NPS Institutional Archive 
□Space Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2017-09 

Photonic compressed sensing Nyquist folding receiver 


Shmel, Richard N. 

Monterey, California: Naval Postgraduate School 


http://hdl.handle.net/10945/56177 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p://w ww. n ps. e du/l ib ra ry 


Caflwuo is the Naval Postgraduate School's public access digital repository for 
research mate rials and institutiional publicatkins created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Caftiouo, NPS's first 
appointed — and published — schoteily author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 Univefsity Circle 
Monterey, California USA 93943 







NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY, CALIFORNIA 

THESIS 


PHOTONIC COMPRESSED SENSING NYQUIST 


FOLDING RECEIVER 


by 


Richard N. Shmel 


September 2017 

Thesis Advisor: 

Phillip E. Pace 

Second Reader: 

Roberto Cristi 


Approved for public release. Distribution is unlimited. 




THIS PAGE INTENTIONALLY LEET BLANK 



REPORT DOCUMENTATION PAGE 

Form Approved 0MB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instruction, 
searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments 
regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden to Washington 
headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and 
to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) Washington DC 20503. 

1. AGENCY USE ONLY (Leave Blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

September 2017 Master’s Thesis 09-25-2016 to 09-15-2017 

4. TITLE AND SUBTITLE 

PHOTONIC COMPRESSED SENSING NYQUIST FOLDING RECEIVER 

5. FUNDING NUMBERS 

6. AUTHOR(S) 

Richard N. Shmel 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 

Monterey, CA 93943 

8. PERFORMING ORGANIZATION REPORT 
NUMBER 

9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

L3 Mission Integration 

10001 Jack Finney Blvd., Greenville, TX 75402 

10. SPONSORING / MONITORING 

AGENCY REPORT NUMBER 

11. SUPPLEMENTARY NOTES 

The views expressed in this document are those of the author and do not reflect the official policy or position of the Department of 
Defense or the U.S. Government. IRB Protocol Number: N/A. 

12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Approved for public release. Distribution is unlimited. 

12b. DISTRIBUTION CODE 

13. ABSTRACT (maximum 200 words) 

The use of integrated optical components and how they can be used to implement a compressed sensing microwave photonic receiver 
are investigated in this thesis. The scope of this research covers the design, simulation, and hardware build of a photonic compressed 
sensing Nyquist folding receiver that is able to directly undersample wideband RE signals and detect the original Nyquist zone and 
frequency information. This is achieved by generating a frequency modulated optical impulse train for directly undersampling the 
RE environment at the antenna. The impulse train undersamples the signals using an optical modulator configuration at 1550 nm 
and collects the detected samples in a low pass interpolation filter. Two independent photonic receiver architectures are designed and 
analyzed over the course of this research. Both receiver designs are simulated in the OPTSIM photonic design suite and are able to 
successfully extract the undersampled signals. The first receiver design proved more effective at representing the signal environment 
and was constructed using electro-optical hardware. Using off-the-shelf components in a proof-of-concept hardware prototype, we 
confirmed the simulation results by correctly identifying input test frequencies up to an order of magnitude higher than the sampling 
frequency. 

14. SUBJECT TERMS 

compressed sensing, undersampling, photonic modulators, NYFR, Nyquist folding, electronic warfare, elec¬ 
tronic support, OPTSIM 

15. NUMBER OF 

PAGES 83 

16. PRICE CODE 

17. SECURITY CLASSIFICATION 

OF REPORT 

Unclassified 

18. SECURITY CLASSIFICATION 

OF THIS PAGE 

Unclassified 

19. SECURITY CLASSIFICATION 

OF ABSTRACT 

Unclassified 

20. LIMITATION OF 
ABSTRACT 

UU 

NSN 7540-01-280-5500 Standard 

Form 298 (Rev. 2-89) 


Prescribed by ANSI Std. 239-18 


1 




























THIS PAGE INTENTIONALLY LEET BLANK 


11 



Approved for public release. Distribution is unlimited. 


PHOTONIC COMPRESSED SENSING NYQUIST FOLDING RECEIVER 


Richard N. Shmel 

Second Lieutenant, United States Army 
B.S., United States Military Academy, 2016 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 

from the 

NAVAL POSTGRADUATE SCHOOL 
September 2017 


Approved by: Phillip E. Pace 

Thesis Advisor 


Roberto Cristi 
Second Reader 


Clark Robertson 

Chair, Department of Electrical and Computer Engineering 



THIS PAGE INTENTIONALLY LEET BLANK 


IV 



ABSTRACT 


The use of integrated optical components and how they can be used to implement a com¬ 
pressed sensing microwave photonic receiver are investigated in this thesis. The scope of 
this research covers the design, simulation, and hardware build of a photonic compressed 
sensing Nyquist folding receiver that is able to directly undersample wideband RF sig¬ 
nals and detect the original Nyquist zone and frequency information. This is achieved by 
generating a frequency modulated optical impulse train for directly undersampling the RF 
environment at the antenna. The impulse train undersamples the signals using an opti¬ 
cal modulator configuration at 1550 nm and collects the detected samples in a low pass 
interpolation filter. Two independent photonic receiver architectures are designed and ana¬ 
lyzed over the course of this research. Both receiver designs are simulated in the OPTSIM 
photonic design suite and are able to successfully extract the undersampled signals. The 
first receiver design proved more effective at representing the signal environment and was 
constructed using electro-optical hardware. Using off-the-shelf components in a proof-of- 
concept hardware prototype, we confirmed the simulation results by correctly identifying 
input test frequencies up to an order of magnitude higher than the sampling frequency. 
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CHAPTER 1: 
Introduction 


Wideband spectrum analysis of electromagnetic waveforms is constrained by many 
factors. Traditional receivers are limited by sampling speeds that must meet the Nyquist 
criterion and components with limited bandwidths, both of which prevent wideband inter¬ 
ception of the radio frequency (RF) environment. The design, simulation, and prototyping 
of a photonic Nyquist folding receiver that alleviates many of the limiting factors that restrict 
current wideband signal collection are discussed in this thesis. 

1.1 Background and Motivation 

Electronic intelligence (ELINT) is intelligence gathering using electronic sensors 
focused on RE signals-of-interest. Target signals vary wildly, ranging from radar, navigation 
beacons, jammers, radio control links, guidance systems, and identification friend or foe 
(lEE) systems [ 1 ]. Signals intelligence (SIGINT) is related to EEINT but is primarily focused 
on communication signals. These terms are often used interchangeably in electronic warfare 
(EW) literature. EW Receivers measure pulsed and low probability of intercept (EPI) signals 
and determine their carrier frequencies, pulse characteristics, phase or frequency variations, 
and polarization [1]. The Nyquist folding receiver (NYER) described in this thesis is able to 
provide pulse rate, carrier frequency, phase variation, and frequency variation information 
while utilizing a compressed sensing methodology. 

The primary motivation for this thesis is to apply compressed sensing techniques to 
a new receiver design implemented with wideband integrated optical components. Progress 
with photonic devices is rapidly surpassing conventional electronic systems. Digital EW 
receivers are taking advantage of these advances in order to digitize and process wide spec¬ 
trum bandwidths. Specifically, optical receivers have an advantage in bandwidth, sampling 
speed, and immunity to electromagnetic interference. Photonic EiNbOs modulators now 
have bandwidths on the order of 40 GHz and can efficiently couple RE-antenna signals di¬ 
rectly into the optical domain [2]. Mode-locked lasers used for sampling now have tunable 
femtosecond-wide pulse widths with pulse-repetition frequencies (PRE) on the order of 300 
GS/s to allow over sampling of complex modulations [3]. 
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1.2 Principal Contribution 

The research objective for this thesis is to determine if integrated photonic compo¬ 
nents can be used to implement a NYFR compressed sensing receiver configured in the 
optical domain as well as characterize the functionality, performance, and limitations of 
these types of photonic receiver architectures. The principal contribution of this thesis is 
the development of a novel compressed sensing receiver that takes advantage of optical 
signal processing techniques. Detailed in this thesis is the entire engineering process from 
the initial architectural designs, through computer simulations, and ending in a hardware 
prototype of the concept. Additional contributions of the research include developing a 
robust method of generating frequency modulated optical pulses and analyzing the func¬ 
tional limitations of photonic sampling. Lastly, the construction of a hardware prototype 
takes the research beyond computer simulations and confirms the physical operation of the 
receiver design. The research methodology was broken up into three primary research 
blocks spanning the designs, simulations, and a hardware prototype build. 

In the first research block, we focus on developing a feasible compressed sensing 
architecture in the optical domain. The NYFR receiver relies on frequency swept optical 
pulses to function as a sampling impulse train. Initially, two separate designs were investi¬ 
gated for implementing a modulated sampling impulse function. The first design examined 
using two separate optical modulators coupled together with one modulator generating an 
optical sampling impulse train and the other incepting the target RF signal. Sampling 
pulses are generated electrically using an electrical pulse generator which is connected to 
a modulator to convert electrical impulses into optical impulses. The sampling frequency 
modulated (FM) Dirac function is coupled to a second modulator that incepts the antenna 
signal onto the optical pulses. This novel sampling method utilizing two modulators is a 
principal thesis contribution. After sampling the RF environment, the optical waveform 
is coupled to a photodiode and amplifier for envelope detection and later digitization. A 
second architecture developed during the research utilizes a single modulator to incept the 
RF signal which is sampled by an analog-to-digital converter (ADC). The trigger for the 
ADC is modulated in order to achieve non-uniform sampling. 

Next, both designs are simulated in the OPTSIM photonic design suite in order to 
confirm the functionality of the NYFR compressed sensing receiver. This was a critical 
step because it allowed the designs to be quickly tested and their limitations confirmed. The 
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test results using simulated RF signals showed that the designs were successful in unfolding 
undersampled waveforms in the optical domain, establishing a foundation for future work 
and confirming the receiver’s effectiveness. Test results indicated that the two modulator 
architecture had superior performance. 

After successful undersampling with RF signals above the Nyquist rate, the next 
block of research focused on a prototype hardware build that could function as a proof- 
of-concept technology demonstrator. Translating a theoretical model into a functional 
receiver required designing a novel electrical pulse-generator circuit that outputs a pulse 
train with a frequency modulated pulse-to-pulse interval. After successfully incepting 
imitation RF signals from a signal generator, the receiver output was sampled and analyzed 
using digital signal processing (DSP) spectrogram techniques. The signals were analyzed in 
a time-frequency representation. Significant effort went into de-noising and edge detecting 
the signals. Hardware tests confirmed the simulation results by correctly unfolding the 
undersampled signals and is a significant research contribution. Ultimately, the hardware 
design proved successful in incepting and extracting RF signals and providing a foundation 
for future work. 


1.3 Thesis Organization 

This thesis report is organized as follows. The NYFR architecture and governing 
equations are described in Chapter 2 along with the technical literature relevant to the thesis 
research. The two photonic NYFR architectures developed for the thesis and the integrated 
photonic components used in designing the receivers are detailed in Chapter 3. Simulation 
models and test results which confirmed the theoretical design functionality are shown in 
Chapter 4. The physical hardware build is described in Chapter 5, along with individual 
component testing. Selected results from the hardware prototype digital signal processing 
steps are highlighted in Chapter 6. Lastly, the thesis research is summarized in Chapter 7 
along with potential areas for future work. The MATLAB code used to unfold the NYFR 
signals and the image processing scripts that were used for Chapter 6 are listed in the 
Appendix. 
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CHAPTER 2: 

Nyquist Folding Methodology 


The model of the photonic receiver presented in this thesis is a compressed sensing 
method known as Nyquist folding. The relevant NYFR theory, governing equations, and 
operation is described in this chapter. In order to effectively sample a signal using only 
one piece of information, the sampling frequency must be greater than twice the maximum 
signal bandwidth. This is the Nyquist criterion and is written as 


fsam ^ ( 2 . 1 ) 

where fsam is the sampling rate, and Bsig is the signal bandwidth. Sampling at a rate lower 
than fsam is called undersampling and causes aliased signals that are symmetrical around 
fsam!'^- The NYFR takes advantage of aliasing in order to fold large bandwidths of the RF 
spectrum into undersampled narrow bandwidths which are later recovered. 

2.1 Analog-to-Information Receiver 

The NYFR functional block diagram is shown in Figure 2.1(a) along with an illus¬ 
tration of the signed Nyquist zones in Figure 2.1(b). The NYFR uses a wideband front-end 
filter to select RF frequencies of interest from the environment. The passband must be 
selected so that the bandwidth is less than or equal to the NYFR maximum intercept fre¬ 
quency. Signals from the wideband filter are sampled by a FM swept clock which is coupled 
to a narrow band interpolation low pass filter (LPF). Lastly, a uniform sampling analog- 
to-digital converter digitizes the signal. Unlike traditional receivers which use a uniform 
sampling rate, the NYFR uses an undersampling rate many times lower than the bandwidth 
of interest [4]. 

Target RF signals are aliased into lower frequency zones which are symmetrical 
around intervals of the undersampling rate cOs- These are referred to as Nyquist zones and 
are diagrammed in Figure 2.1(b). The Nyquist zones are named in an alternating sign 
convention and are each (jOsI2 rad/s wide. If the sampling rate is uniform, these aliased 
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signals are indistinguishable from the true center frequency; however, the NYFR samples 
with a modulated linear chirp slope. The sampling waveform imparts a modulation on the 
target signal which is dependent on the original Nyquist zone of origin [5]. 


(a) 



(b) 

Signed Signal Nyquist Zone N 



Analog Frequency 

Figure 2.1. Analog-to-information Receiver Showing a) NYFR Operation and 
b) Signed Nyquist Zones. Adapted from [4]. 


The FM sampling clock generates a sampling signal centered around an undersampling 
frequency cOs with a RF sampling clock modulation 6(t) as 


(fit) = cOst + 6{t) 


( 2 . 2 ) 


and is used to sample the RF spectrum. An example target signal is modeled as 
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x{t) = cos [oJct + <A(0) 


(2.3) 


where cOc is the eenter frequeney and ij/ represents the RF signal information eontent. 
Convolving the two signals through the LPF yields the interpolation filter output 


y(t) = cos {a>ft +- M6(t)) (2.4) 

where a>f is the output folded frequeney in Nyquist zone N = 0, /3 = sgn{a)c - cosk) and 
is the folded speetral orientation, and M = f3k whieh is the modulation seale faetor [5]. 
Also, k = round[tOclcOs) and is ealled the sampling harmonie faetor [4], [5]. The indueed 
modulation seale faetor M is measurable and 6(t) is known, so the original signal Nyquist 
zone ean be estimated. Expanding tof from (2.4) results in the output folded frequeney 


tOf = \tOc- tOsk\ (2.5) 

whieh shows the relationship between the sampling modulation and RF eenter frequeney [5]. 

2.2 Unfolding the Undersampled Signals 

In order to extraet original frequeney information from the aliased signals in Nyquist 
zone N = 0, they must be unfolded and the ambiguities resolved. From (2.4) and (2.5) we 
see that the modulated RF signal frequeney is a funetion of the original eenter frequeney 
and the applied ehirp modulation. The minimum differenee between the signal modulation 
and eloek modulation sealed by the Nyquist zone indieates the signal’s original frequeney 
region. The target frequeney is extraeted from the aliased signal in the N = 0 Nyquist zone 
by first generating an array Z[i] of the modulated signal slope Sm minus the FM eloek slope 
Sc sealed by the signed Nyquist zones M(i). The array is modeled as 

Z[i] = {S^ - Mii)Scf (2.6) 

where Sc is the slope of the eloek modulation, Sm is the slope of the pulse train in time- 
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frequency analysis, and M{i) is the signed Nyquist zones N = [0, -1,2, -3...], respectively. 
The next step is to find the index i of the minimum value of the array so that i = min[Z\i\). 
Lastly, the unfolded Nyquist zone is defined by 


N = sgn{M{i)){i - 1) (2.7) 

where N is the original signed Nyquist zone of the RF signal. Adding the center frequency 
from Nyquist zone N = 0 yields coc [6]. This is shown in Section 4.1 for a 2.75-GHz test 
signal applied to the photonic NYFR simulation. 

The NYFR compressed sensing theory and operation was detailed in this chapter 
by first modeling the NYFR architecture and then describing the governing equations. The 
analog-to-information process uses clock modulation to impart a slope onto an undersam¬ 
pled signal, which is then used to unfold the aliased waveform. The photonic NYFR designs 
as well as the integrated optical components used in the architecture are detailed in the next 
chapter. 
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CHAPTER 3: 

Receiver Operation and Design 


The problem presented in this research is the implementation of a photonic NYFR 
prototype. Correct implementation first requires a detailed design and analysis of required 
components. Additionally, operating limits and constraints must be taken into account 
during the design phase. In order to test the feasibility of a photonic compressed sensing 
receiver, two independent designs are investigated. The first is a purely optical design 
intended as the primary research contribution of this thesis. The second receiver is an 
electro-optical solution that functions as an interim design and technology demonstrator for 
the proposed architecture. 

The design phase of the research and the integrated optical components used in the 
photonic NYFR are described in this chapter. The two receivers are discussed in detail 
with emphasis on design choices and operating methodology. The constraints of optical 
sampling components along with each receiver’s operating limitations are investigated in 
the last section. 


3.1 Integrated Optical Components 

The optical components used to design the NYFR are described in this section. Both 
design solutions integrate laser sources, optical modulators, and optical detectors in order 
to undersample and modulate a signal waveform. The specific component types and their 
operations are expanded upon in subsections. 

3.1.1 Semiconductor Laser Source 

A continuous wave (CW) semiconductor laser is used as the optical emitter for the 
NYFR. The specific design is a buried heterostructure laser consisting of a gallium arsenide 
(GaAs) bandgap junction embedded in an AlGaAs substrate. A wavelength of 1550 nm 
was used to limit attenuation over the connecting single mode fiber optic cable. The laser 
operates by applying current above a threshold value 4, increasing the probability of carrier 
recombination in the GaAs bandgap and resulting in stimulated emission. The emitted 
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optical power is written as 


Pe 


- nv 

e 


(3.1) 


where I is the applied eleetrie eurrent, h is Planck’s constant, v is the frequeney, e is 
elementary eharge, and rji is the probability of earrier reeombination in the gap region [7]. 
In cases where the laser is modulated direetly through the supply current, the applied eleetrie 
current has both AC and DC eomponents [7] 


I = Io+ Zoe'"'”' (3.2) 

where cOm is the modulation frequeney and Iq and z’o are eonstant DC and AC eurrent 
magnitudes. 

3.1.2 Optical Modulators 

Both reeeiver designs used one or more eleetro-optie modulators to modulate a 
CW laser source. Mach-Zehnder Interferometers (MZIs) are used to facilitate amplitude 
modulation through splitting and reeombination of in-phase and out-of-phase beams. The 
photonie NYFR designs used single-arm MZIs eomposed of a metallie-infused lithium 
niobate waveguide ehannel which splits into two branches at a Y-junction. The upper braneh 
is braeketed by eleetrodes whieh apply a voltage potential to the upper waveguide. Voltage 
applied to the eleetrodes creates an eleetromagnetie field through the niobate strueture. This 
alters the propagation constant on the braneh and ehanges the traveling wave’s phase. A 
seeond Y-junetion reeombines the two branehes. Depending on the phase difference A0, 
we see that the two waves undergo either eonstruetive or destruetive interference whieh 
modulates the output optical intensity [8]. This is illustrated in Figures 3.1(a) and 3.1(b) 
for an applied eleetrode voltage 14. 
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Figure 3.1. Diagram of Single-Arm MZI Showing a) In-Phase 
Constructive Interference and b) Out-of-Phase Destructive Interference. 

Adapted from [8]. 

The output optical intensity Ig of an MZI is a function of the A0 path length phase 
difference and is 


lo 


h 

2 


I -I- cos (A<p) 


(3.3) 


where 7/ is the incident optical intensity. The phase difference between the recombining 
waveforms from each MZI branch is 


II 
















































(3.4) 


A. / 

A0 = \ n— + (po 

where V^r is the switehing/threshold voltage, Va is the applied eleetrode voltage, and ipo is 
an initial phase offset that is built into the device or created with a DC voltage known as a 
DC bias [9]. When there is a phase difference of A0 = n, the output has its lowest value, 
while a phase difference of zero results in the highest output value. Combining (3.3) and 
(3.4), we obtain the expanded MZI output intensity 


lo 


h_ 

2 


1 + cos 


n -h 



(3.5) 


The modulator input/output intensity relationship, or transfer function, is represented 
graphically in Figure 3.2. 



Electrode voltage (vj 


Figure 3.2. Graphic Representation of MZI Function with Both the Quadra¬ 
ture and DC Bias Marked 
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For a modulator with equal guided-wave path lengths, the application of a DC bias translates 
the curve left and right. The use of an attenuator or amplifier at the antenna increases or 
decreases the value of 14, respectively. This shows the output intensity as a function of the 
electrode voltage and illustrates the 14 switching voltage. The quadrature point is created 
by an initial phase constant oicpo = n manufactured into the MZI so that any voltage applied 
shifts the output in the most linear region of operation [10]. This configuration is ideal 
for coupling an AC small signal onto an optical carrier. The effect of a DC bias is shown 
as well, where a DC voltage is used to shift the curve such that the initial 14 = 0 point is 
perfectly destructive interference and V = I 4 is the peak intensity point. This configuration 
is used for optical pulse generation. 

3.1.3 Optical Detectors 

The photonic NYFR designs use PIN photodiode detectors to convert optical signals 
into electrical signals. These detectors consist of p-type and n-type doped silicon separated 
by an intrinsic region of high resistivity material. Absorption of photons in the silicon 
material causes electrons to drift and creates current [11]. The NYFR photodiodes are 
matched to the 1550-nm wavelength of the CW laser. The two most important factors for 
PIN diodes are gain and frequency response. The NYFR photodiodes have an internal low 
noise amplifier (LNA) with a 10 dB gain and a frequency response bandwidth of 20 GHz. 
This response is primarily limited by the transit time of charge carriers separating from the 
substrate and moving to opposite ends of the depletion region. An internal shunting effect 
is also present and is calculated as [11] 


OJp = 


1 


(3.6) 


where ojp is the photodiode cutoff frequency in rad/s and Ri and Q are the internal resistance 
and capacitance of the diode, respectively. 


3.2 Double-Modulator Design Concept 

The double-modulator NYFR (DM-NYFR) architecture is an implementation of a 
photonic compressed sensing receiver. It is designed around two MZI modulators. The first 
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functions as a sampling impulse train generator while the second incepts the RF bandwidth. 
In the photonic NYFR architecture, the CW laser functions as a carrier for the signal and is 
analogous to electric current on a wire. 

3.2.1 DM-NYFR Operation and Block Diagram 

The overall block diagram of the DM-NYFR is shown in Figure 3.3. The front end 
of the architecture is a fiber-coupled CW laser at 1550 nm. This is conneeted to the first MZI 
whieh converts the CW laser into a sampling pulse train with a linearly swept pulse repetition 
frequency (PRF). An electrical signal P{t) is connected to the single-arm eleetrodes of MZI 1 
and drives the modulation scheme by linearly sweeping a short electrical impulse train. The 
output optieal pulse intensity, represented as 4, is a function of the frequency and amplitude 
of the electric signal as well as the MZI transmission (3.5). The front end optical pulse 
generation is discussed in detail in Section 3.2.2. 



Figure 3.3. DM-NYFR Block Diagram Showing the Overall Architecture of 
the Receiver 


The FM pulse train functions as a sampling signal. When coupled into the second 
MZI with the RF-signal input, the short-time duration pulses emulate the sampling eleetrical 
signal from (2.2). As the FM optieal pulse train is passed into MZI2, RF signals from the 
environment are amplified from an antenna and used to drive MZI2’s electrode inputs. The 
amplitude of the input optical pulses from MZIl are modulated with the amplified antenna 
signal voltage, encoding the RF signals as an amplitude modulated (AM) envelope on the 
sampling pulses. This eoneept is illustrated in Figure 3.4, in which a diagram of an arbitrary 
sinusoidal signal encoded on a sampling pulse train is shown. This is a good representation 
of MZI2’s optical output. 
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Figure 3.4. Example of an Arbitrary Signal Sampled by an Optical Pulse 
Train 


The output from MZI2 is coupled into a PIN photodiode with an internal LNA in 
order to convert the optical signal into electrical current. An interpolation low pass filter 
(LPF) is connected in series to the detector output. The output of the interpolation filter from 
the photonic NYFR is equivalent to y(t) (2.4). The LPF bandwidth must be greater than 
or equal to cosll, which is the first Nyquist zone upper frequency shown in Figure 2.1(b). 
A sampling ADC driven by a clock input at IcOs digitizes the information at the required 
Nyquist rate. Lastly, the ADC passes this information to the digital signal processing (DSP) 
block where the data is de-noised and unfolded to extract the RF information from the 
environment. 

3.2.2 Optical Pulse Generation 

The front end optical pulse generation using MZIl shown in Figure 3.3 is the most 
important part of the DM-NYFR architecture. Generating FM optical pulses with small 
pulse widths enables the NYFR to differentially alias signals to lower Nyquist zones for 
detection. In the DM-NYFR, optical pulses are generated by applying a pulse signal P{t) 
to MZIl so that the coupled CW laser is selectively switched at a linear chirp rate. The 
electrical pulse signal at the input is modeled as an FM Dirac comb impulse train, comprised 
of n number of 6{t) delta functions. This is an approximation, because the Dirac comb 
impulses have zero pulse width, which is impossible to achieve in actual operation. If the 
pulse width is much smaller then the RF signal’s period, a zero-width impulse train is a 
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good approximation. As a result, P{t) is represented as a sum of impulses 

P{t) = V,Yj^{t-n{t)) (3.7) 

where 14 is the voltage amplitude of the MZI switching voltage and n is the impulse trigger. 
In order to model the variable PRF of P{t), the delta functions must be triggered from a 
linear FM sinusoid. In the DM-NYFR architecture, this is modeled as a chirp signal Sc{t), 
with n{t) set on every positive-to-negative zero crossing such that 

n{t) = ZCR~{Sc{t)) (3.8) 

where ZCR~ is the notation for a single-ended zero crossing. For the DM-NYFR, two types 
of linear chirp signals are used, a sawtooth waveform Sci(0 and a triangular waveform 
Sciit)- The sawtooth chirp is given by 


Sci{t) = sin2;r|/ci + (3.9) 

where fc is the center frequency, AF is the chirp bandwidth, and tm is the chirp period. The 
triangular waveform is defined as 


5 ^2(0 = sin27r[5i(0 + ^2(0] (3.10) 

with Si as the rising chirp sweep and S 2 as the falling chirp sweep. The chirp sweep 
functions are given by 


Si(o = 



AF , 


(3.11) 


for the positive slope and 
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(3.12) 


S2it) = 



AF 



2 


for the negative slope, where fc is the center frequency, AF is the chirp bandwidth, and is 
the chirp period [12]. The linear chirp sinusoid (3.9) plotted along with the pulse function 
(3.8) is detailed in Figure 3.5. Setting the amplitude of the pulse signal to K- switches 
the optical modulator at the pulse repetition frequency (PRF) defined by (3.7) and coverts 
the input CW laser into a sampling impulse train. This illustrates how the front-end pulse 
generation block works in the DM-NYFR architecture. 



T 


Figure 3.5. An Example Linear Chirp Used to Generate a FM Dirac Comb 
Every Positive-to-Negative Zero Crossing 


3.3 Single-Modulator Design Concept 

The single modulator design that serves as a simplified technology demonstrator for 
the photonic NYFR is described in this section. Conceived as an alternative architecture 
early in the research process the single-MZI NYFR (SM-NYFR) is an electro-optical hybrid 
solution that takes advantage of the high bandwidth of photonic components but with the 
low cost and complexity of electronic sampling devices. The trade off in the electro-optical 
design is a lower maximum RF-interception frequency than the DM-NYFR solution and 
higher receiver noise level. The SM-NYFR block diagram is shown in Figure 3.6. 
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The SM-NYFR uses a CW laser eoupled to a single-arm MZI biased at quadrature. 
The eleetrodes are attaehed to an antenna/LNA pair whieh ineepts RF signals from the 
environment. The high bandwidth of optical components, greater than 40 GHz, ensures 
that a large RF spectrum is modulated onto the laser. The MZI output is connected to a 
PIN photodiode, which converts the photonic signal into an electrical signal. The signal is 
filtered with an interpolating LPF that has a bandwidth greater than the maximum frequency 
of the first Nyquist zone. An ADC samples the electrical signal driven by a clock pulse P{t) 
detailed in (3.7). The digital information is processed by a DSP, which unfolds and extracts 
the RF frequencies from the environment. 




Antenna Input 


Legend: 

- Optical signal 



Figure 3.6. SM-NYFR Block Diagram Showing the Overall Architecture of 
the Receiver 


3.4 Design Constraints and Considerations 

The photonic NYFR’s primary operating boundaries and constraints come from 
the optical sampling with interferometers. The optical pulse width and ADC resolution 
determine the maximum frequency that can be modulated onto a laser source. Both pulse 
width and sampling limitations are examined in detail in this section. The performance of 
potential hardware builds is also analyzed in order to determine the maximum RF-signal 
intercept frequencies. 

3.4.1 Pulse Width Limitations 

The photonic NYFR relies on signal sampling by optical pulses where the instan¬ 
taneous power of the pulse train is proportional to the sampled waveform. In order to 
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effectively sample a target RF signal, the pulse width must be smaller than the temporal 
variation of the modulating signal [13]. The 3 dB bandwidth Bp of an optical pulse is 
the inverse of the pulse width Tp or Bp = IjTp. From the Nyquist criterion, the sampling 
frequency fs for an arbitrary signal with bandwidth B is fs = 2B, written in terms of pulse 
width as 


< I (3.13) 

where Tp is the largest allowable pulse width that can effectively sample a signal with a 
period of Ts. For the photonic NYFR design, this means that optical pulse widths limit the 
maximum frequency that can be sampled by MZI2, and pulses must be smaller than half 
the largest RF signal period. 

3.4.2 Electro-Optical Sampling Limitations 

The optically sampled signal from MZI2 is coupled into a PIN photodiode which 
converts optical power into an electrical signal. The output voltage must then be quantized by 
a sampling ADC with finite resolution. Digitizing an analog waveform imparts a sampling 
error which constrains the photonic NYFR’s operating limits. Let the signal applied to the 
MZI electrodes be 


V{t) = Asin(2;r/mt) (3.14) 

where fm is the maximum signal frequency and A is the signal amplitude [14]. For a 
sampling pulse width Tp, the error in voltage 6V for an arbitrary signal at time in a 
traveling wave optical device is represented by 

p_Tp 

6V=^f A sm(2n - A sinilnfmti). (3.15) 

Tp 4 +^ 

The voltage error is the difference in instantaneous voltage integrated over the width of 
a sampling pulse. Using a Taylor series expansion to evaluate the integral (3.15) for the 
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maximum voltage error yields 


(nfmTpfA 

\SVU. = / • (3.16) 

o 

For a n-bit sampling ADC, the voltage error must be less than half of the level spaeing 
whieh is 


\^V\inax ^ 


2A n \ _ A 

2"\2I ~ 2"' 


(3.17) 


This results in a Tp pulse width limitation of 


Tp< 


2«-i 


1/2 


^ fm 


(3.18) 


where fm is the maximum target frequeney and n is the resolution of the ADC [14]. 


3.4.3 Implications for Photonic NYFR Design 

The eleetro-optieal sampling error assoeiated with the ADC is the dominate limita¬ 
tion for the photonie NYFR arehitecture. Design implieations are shown in Figures 3.7 and 
3.8, whieh are plots of the maximum intereept frequeney as a funetion of pulse width for 
six and eight bit ADCs, as well as for two separate hardware designs. Figure 3.7 is an illus¬ 
tration of a prototype build using inexpensive off-the-shelf eomponents with pulse widths 
in the nanoseeond range. The upper limit build illustrating pieoseeond and sub-pieoseeond 
pulse widths, whieh are aehievable with either high quality eleetrieal pulse generators or 
mode loeked fiber lasers, is shown in Figure 3.8. There is a elear design trade off between 
maximum intereept frequeney and ADC resolution, with higher resolutions eorresponding 
to lower target RF signals. Additionally, higher frequeneies require shorter pulse widths, 
whieh inereases eosts. 
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Figure 3.7. Maximum Intercept Frequency for Nanosecond Pulse-Widths 
Illustrating a Low-Cost Prototype Build 



Pulse width (ps) 


Figure 3.8. Maximum Intercept Frequency for Picosecond and Sub- 
Picosecond Pulse-Widths Illustrating a Potential Fligh-Cost Build 


In summary, both photonic NYFR designs and their operating methodologies were 
diseussed in this ehapter. The eonstraints, limitations, and design trade ofFs for the photonie 
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NYFR architecture are examined, and modeling shows that the optical sampling error and 
the resulting optical pulse width is the primary limiting factor. Additionally, the integrated 
optical components used in the NYFR are were discussed. The computer simulations and 
results for both photonic NYFR designs are covered in the next chapter. 
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CHAPTER 4: 

Photonic NYFR Computer Simulations 


The computer simulations of the two receiver designs are described in this chapter. 
Both the DM-NYFR and the SM-NYFR architectures are modeled and simulated in order 
to test their operation. The final models for both the double and single modulator designs, 
as well the key test results for each architecture, are detailed in this chapter. Simulation 
test results were successful in demonstrating correct inception and unfolding of target RF 
signals. 


4.1 DM-NYFR Computer Simulation 

The DM-NYFR computer model is shown in Figure 4.1 while Table 4.1, lists the 
model components and their functions. 


Antenna input 



Figure 4.1. DM-NYFR Computer Simulation Model Showing All Compo¬ 
nents and Connections. Specific Functional Sections are Flighlighted and 
Labeled with Dashed Boxes. 
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Table 4.1. Component Labels, Names, and Descriptions from the 
DM-NYFR Computer Simulation Model 


Label 

Name 

Function 

SAWTOOl 

Sawtooth wave generator 

Generates a electrical sawtooth wave from Vmin to 

Vmax that is fed into the VCO. 

VCOl 

Voltage controlled oscillator 

Generates a linearly chirped sinusoid from the saw¬ 
tooth input based on voltage sweep. 

ZCGENl 

Zero crossing generator 

Generates a picosecond wide electrical pulse every 
positive-to-negative zero crossing. 

egain2 

Electrical gain 

Amplify IV pulse to 14 (5.5V). 

scope 1 

Electrical scope 

Electrical oscilloscope measuring tool. 

scope2 

Electrical scope 

Electrical oscilloscope measuring tool. 

laslorl 

Easer source 

1550nm CW laser source. 

loamodl 

MZIl 

Eirst optical modulator, uses EM electrical pulses to 
generate EM optical pulses. 

probe1 

Optical probe 

Optical signal measuring tool. 

photod_pinl 

PIN Photodiode 

Converts an optical signal into an electrical signal. 

egainl 

Electrical gain 

Amplify electrical signal from the photodiode. 

scope3 

Electrical scope 

Electrical oscilloscope measuring tool. 

SIGGENl 

Signal generator 

Generates an electrical target signal for the NYER 

SIGGEN2 

Signal generator 

Generates an electrical target signal for the NYER 

UNIERR2 

Uniform noise 

Adds noise to the signal input 

loamod2 

MZI2 

Modulates the EM pulse train to sample the antenna 
signal. 

probe1 

Optical probe 

Optical signal measuring tool. 

photod_pin2 

PIN Photodiode 

Converts an optical signal into an electrical signal. 

egainS 

Electrical gain 

Amplify electrical signal from the photodiode. 

filbutl 

Butter worth EPE 

8th order butterworth EPE set to 2/^. 

scoped 

Electrical scope 

Electrical oscilloscope measuring tool. 


The DM-NYFR design uses two coupled MZIs to both sample and incept a RF signal 
from the environment. The most important aspect of the design is a front-end optical pulse 
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generation block that converts a CW laser into a FM optical pulse train. The computer model 
is broken up into functional sections that each perform a specific task. As shown in Figure 
4.1, the sections are electrical pulse generation, optical pulse generation, signal modulation, 
antenna input, and sampling. The rSoft OPTSIM design suite from Synopsys is used as a 
photonics modeling environment. This software package enables high fidelity simulation of 
integrated optical systems, and both receiver designs were successfully simulated and tested. 
OPTSIM uses a block simulation design environment where components and functions are 
modeled as discrete units that are linked in a user interface. 

4.1.1 Electrical Pulse Generation 

The electrical pulse generation block models (3.8) in the rSoft design suite and 
functions as a proof-of-concept for a hardware build. An electrical signal generator outputs 
a sawtooth waveform, varying the amplitude from Vmin to Vmax over a modulation period 
tm- The signal is sent to a voltage controlled oscillator (VCO) that chirps a sinusoid wave 
from fmin to fmax- For the VCO operation, V^in maps to and V^ax maps to f^ax- The 
modulation period matches tm- The chirped sinusoidal wave from the VCO is sent to a zero 
crossing pulse generator that outputs a 1 .OV electrical pulse with a picosecond-wide pulse 
width every positive-to-negative transition. In order to modulate the CW laser, the l.OV 
pulse train is amplified to the MZI switching voltage 14- This is shown in Figure 4.2. 

(a) (b) 




Figure 4.2. Electrical Pulse Generation Simulation Outputs With A) Input 
Sawtooth Waveform and B) Output Pulse Train Frequency Representation 
Showing FM Chirp 


25 







































Figure 4.2(a) is an illustration of the input sawtooth waveform as a time-amplitude repre¬ 
sentation from seope one. Figure 4.2(b) is an illustration of the output bandwidth of the 
pulses sampled at scope two. The x-axis is frequency, and the y-axis is relative power in 
dB. This shows a wide bandwidth AF which is the chirp optical pulse train. 

4.1.2 Optical Pulse Generation 

Electrical pulses are used to drive the first MZI and generate optical pulses as in 
(3.7). The MZI is simulated as a single-arm cosine squared modulator. The FM optical 
pulse train is shown in Figure 4.3, which is a MATLAB spectrogram of the data from probe 
one. Spectrograms use a short time Fourier transform (STFT) window scanned over a signal 
to make a time-frequency representation of the waveform. The optical pulses are chirped 
in a sawtooth pattern from the driver signal on the MZI. The input pulses are centered at 
1.0 GHz and modulated with tm = 5 /us and AF = 150 MHz. The output optical pulses in 
Figure 4.3 match the MZI inputs, indicating a successful modulation technique. 



Time (MS) 

Figure 4.3. Optical Pulse Train Output from MZIl. The Spectrogram 
Shows a Successful FM Pulse Sweep Taken from Probe One. 


4.1.3 Antenna Input and Signal Modulation 

The antenna input is comprised of an adjustable number of sinusoidal signal com¬ 
ponents combined in parallel with a uniform noise generation block. This allows variations 
in the signal-to-noise ratio (SNR) for testing. The input signals are single frequency tones 
due to software limitations. The antenna block drives MZI2, and the input optical pulses 
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sample the RF signals. The resulting signal information is converted from the optical do¬ 
main into the electrical domain and then filtered and collected in an interpolation LPF. The 
LPF output is sampled by scope six with sampling rate 2a)s and piped into MATLAB for 
signal processing. 

4.1.4 DM-NYFR Simulation Results 

The simulated DM-NYFR architecture is tested with target signals in order to see if 
the model can successfully undersample and modulate signals for recovery by the DSP. The 
simulation is conducted with a sawtooth FM modulation. The modulation period tm = 5.0 
fxs with fmin = 1 GHz and fmax = 1-15 GHz resulting in a AF sweep bandwidth of 150 
MHz. From Figure 2.1 this gives Nyquist zones every 500 MHz, with zone N = 0 from 0 
to 500 MHz. The test results show the first two Nyquist zones (A = 0, -1) as well as the 
modulation signal. Scope six, the sampling ADC, has a 2fs sampling rate of 2.0 GHz, and 
the LPF cutoff frequency is 2/^. For all conducted tests, the SNR was 10 dB [15]. 

Test One: No Modulation 

The first test checks the architectures response with no FM sweep applied to the first 
MZI and the results are shown in Figure 4.4. 
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Figure 4.4. DM-NYFR Test One: No Modulation. Disconnected VCO 
Re-sults in Aliased Signals which Are Unresolvable to Their Original 
Frequencies. 


The VCO was disconnected, so the sampling pulse train had a constant frequency of 1.0 
GHz. The applied signal is 2.7 GHz, so the expected results when examining the DSP 
output are aliased signals centered around 500.0 MHz. The output shown in Figure 4.4 
confirms the architecture with no modulation. Aliased signals at 0.3 GHz and 0.7 GHz 
are seen along with the 1.0 GHz sampling signal. Without modulation the undersampled 
signals cannot be resolved to their original frequencies. 

Test Two: 2.75 GHz Signal 

The second test is the architectures response with the FM sweep enabled. The 
optical pulse train has an FM chirp and successfully undersamples the signal. The target 
signal is 2.75 GHz, chosen because that is a common aerial search radar frequency. The 
received signal is shown in Figure 4.5 in the N = 0 Nyquist zone with a modulated slope 
from the MZIl driver signal. The results are a time-frequency spectrogram computed in 
MATLAB. Two Nyquist zones, V = 0, -1, are shown completely, and are mirror images of 
each other. Also visible is sampling signal from 1.0 to 1.15 GHz. 
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Figure 4.5. DM-NYFR Test Two: 2.75 GFIz Signal. Sampling FM Pulses 
Successfully Modulating a Signal Using NYFR Architecture. 


From (2.6) and (2.7), the original frequency can be extracted from the sampling 
signal and the folded signal. A MATLAB program is used to unfold the signal, and 
the results are displayed in Figure 4.6. The lowest difference is Nyquist zone N = -5, 
which corresponds to signals between 2.5 and 3.0 GHz. Estimating a center frequency of 
approximately 0.15 GHz from Figure 4.5 and adding it to the lower frequency bound of the 
Nyquist zone yeilds the correct original frequency of 2.75 GHz. 
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Figure 4.6. Unfolding the Original 2.75 GHz Signal. Using (2.6) to Compute 
the Scaled Lowest Slope Difference. Highlighted Value is the Likely Original 
Frequency Zone. 


Test Three: 5.1 and 5.65 GHz Signals 

The next test was a 5.1 GHz signal and a 5.65 GHz signal, which are some typical 
C-band radar frequencies. Each signal was applied separately and the results analyzed 
in MATLAB. The results are shown in Figure 4.7, which displays both signals in the 
N = 0 Nyquist zone. Comparing the signal slopes and modulation signals, we can extract 
the original frequency ranges. Adding the center frequencies from the results yields the 
original frequencies of the signals. Visible in Figure 4.7 are weaker images that appear to be 
signals; however, these are actually aliased signals from higher Nyquist zones folding down 
due to the wide bandwidth of the interpolation FPF. Reducing the FPF cutoff frequency to 
the upper limit of Nyquist zone N = 0 removes these artifacts. 
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Figure 4.7. Test Three Results Displaying A) 5.10 GHz and B) 5.65 GHz 
Signals. Both Target Waveforms are Successfully Undersampled and Recov¬ 
erable from the Spectrogram. 


Test Four: 2.7 and 9.1 GHz Combination Signal 

This test was conducted to determine how multiple target signals are resolved by the 
DM-NYFR simulation. Two target waveforms were applied to the antenna input: 2.7 GHz 
at 100 mVpp and 9.1 GHz at 50 mV^^. The uniform noise block was disconnected for this 
test. A spectrogram of the instantaneous optical power at MZI2’s output is shown in Figure 
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4.8. Visible is the FM signal at 1.0 GHz and the modulated target signals in Nyquist zone 
N = 0. The antenna signals exhibit non-uniform aliasing expected from the NYFR, and 
two distinct signals are clearly identifiable. Given the slopes and center frequencies from 
the spectrogram, (2.6) and (2.7) can be used to successfully unfold the input RF signals to 
their original values. 



1 2 3 4 5 6 7 8 


Time (ps) 

Figure 4.8. Test Four Results Displaying the Spectrogram of Folded 2.7 GFIz 
and 9.1 GFIz Signals 


4.2 SM-NYFR Computer Simulation 

The SM-NYFR uses a single MZI that acts as a carrier for the target signal. The 
compressed sensing is performed by an ADC triggered with a FM chirp signal described 
by (3.8). The rSoft simulation model is shown in Figure 4.9 and is comprised of an antenna 
block, signal modulation block, optical-to-electrical conversion stage, compressed sensing 
block, and FM sampling pulse generation. Individual components are listed in Table 4.2. 
The FM timing signal for the ADC results in a modulated sampling scheme which allows 
the compressed signal to be unfolded. 
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FM sampling pulse generation 



Figure 4.9. SM-NYFR Simulation Model with all Blocks Outlined 


Table 4.2. SM-NYFR Component Labels, Names, and Descriptions from the 
SM-NYFR Computer Simulation Model. 


Label 

Name 

Function 

SAWTOOl 

Sawtooth wave generator 

Generates a electrical sawtooth wave from Vmin to 

Vmax that is fed into the VCO. 

VCOl 

Voltage controlled oscillator 

Generates a linearly chirped sinusoid from the saw¬ 
tooth input based on voltage sweep. 

ZCGENl 

Zero crossing generator 

Generates a picosecond wide electrical pulse every 
positive-to-negative zero crossing. 

SIGNALl 

Electrical signal 

Target RE signal. 

NOISE 1 

Noise generator 

Generates uniform electrical noise. 

MZI_MODl 

MZIl 

Modulates antenna signal onto CW laser 

EASERl 

Easer source 

1550nm CW laser source. 

PIN_PH0T01 

Photodiode 

Converts an optical signal into an electrical signal. 

SAMPEEl 

ADC sampler 

Converts analog signals to digital signals and is 
driven by ZCGEN1. 
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4.2.1 SM-NYFR Simulation Results 

The simulations were eondueted as deseribed in Seetion 4.1.4. Testing was done to 
determine the performanee of the SM-NYFR in comparison to the DM-NYFR. From test 
results, we see that the SM-NYFR has similar performance but a higher noise level with 
more artifacts in the results. 

Test One: 5.25 GHz RF Signal 

A 5.25 GHz signal was applied to the MZI electrode input and the results captured 
and exported to MATLAB. Two tests are run with the 5.25 GHz signal, the first with the 
VCO disconnected and no modulation present and the second with the VCO outputting a 
FM chirp to the zero-crossing generator to create an FM sampling pulse train used as an 
ADC trigger. The spectrogram of the 5.25 GHz signals are shown in Figure 4.10. From the 
results, we see Nyquist zone N = 0 and the modulated signal can be successfully unfolded 
to the correct frequency. The noise level appears higher than in the SM-NYFR results. 


(a) (b) 



123456789 123456789 

Time (ps) Time (ps) 

Figure 4.10. 5.25 GHz RF Results from an A) Unmodulated Sampling Signal 
and B) Modulated Sampling Signal 


Test Two: 8.10 GHz RF Signal 

For a second test, an 8.10 GHz signal was applied to the antenna input. The result 
is shown in Figure 4.11, and the original frequency was extracted. Unlike the DM-NYFR, 
from results we see large amount of noise clutter that could obscure a low signal-to-noise 
ratio signal such as an LPI waveform. This is indicated in Figure 4.11 by a label which 
marks an ambiguous area where the signal blends into the noise background. 
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Figure 4.11. 8.10 GHz Target Signal Response. The High Noise of the 
SM-NYFR Design Makes Low SNR Signal Detection DifPicult. 


In summary, the computer simulations of the two receiver designs were discussed 
in this chapter. Both the DM-NYFR and the SM-NYFR architectures were modeled 
and simulated in order to test their operation. While both designs correctly unfolded 
compressed, incepted RF signals, the DM-NYFR design has a better noise response. The 
hardware prototype design of the DM-NYFR following the successful simulations is shown 
in the next chapter. 
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CHAPTER 5: 

Receiver Hardware Prototype Build 


After successfully confirming the photonic NYFR designs in rSoft’s OPTSIM de¬ 
sign suite, a prototype NYFR was built to demonstrate physical operation. The design, 
construction, and testing of the hardware NYFR prototype and components are discussed 
in this chapter. Time constraints prevented building both the DM and SM architectures, 
and because the main focus of the research was on implementing integrated photonic com¬ 
pressed sensing techniques, only the DM-NYFR prototype was built in hardware. The 
principal challenge in constructing the hardware prototype was integrating commercial of- 
the-shelf (COTS) components in unique configurations that were not considered when the 
equipment was being manufactured. The hardware build makes use of components never 
intended for SIGINT or FLINT applications. 

Additionally, component cost and funding was also a deciding factor during the 
prototype build. Due to funding limitations, lower cost equipment was used that did not 
meet the intended performance specifications of the designed or simulated architecture. 
Consequently, the prototype build is a scaled-down design that functions as a technology 
demonstrator and proof-of-concept rather than a fieldable receiver. Specifically, although 
the NYFR was designed and simulated to function in the GHz range with Nyquist zone band- 
widths of 500 MHz, the hardware prototype is only able to operate in the kHz/MHz range 
with zone sizes of 50 kHz. Despite these constraints, the prototype is able to successfully 
undersample at sub-Nyquist rates and extract the original frequency information. 

5.1 Overall Setup and Equipment Used 

The DM-NYFR prototype build matches the block diagram in Figure 3.3 and is 
shown in Figure 5.1. The components list is shown in Table 5.1. The prototype uses a 
ThorLabs S1FC1550 1550 nm laser with an output power of 1.75 mW coupled into a JDSU 
X5 single-arm MZI with a DC bias to move the MZI out of positive quadrature as shown in 
Figure 3.2. The DC bias is set so that an applied electrode voltage of O.OV will combine the 
MZI arms out-of-phase, and an electrode voltage of 5.5V will combine the beams in-phase. 
The X5 electrodes are attached to an electrical pulse generator (EPG) driven by an Agilent 
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E3612A DC power supply and Agilent 33220A 20 MHz funetion generator. The EPG is 
detailed in Seetion 5.3. 

The MZIl output is eoupled to the MZI2 input, whieh is also a X5 modulator driven 
with another 33220A funetion generator whieh imitates an antenna input. The MZI2 DC 
bias is not used as the phase eonstant cpo set during manufaeture biases the modulator at 
positive quadrature. This allows the applied AC small signal to aet on the most linear region 
of the MZI transfer function. The MZI2’s output is connected to a Newport D-30IR detector. 
The detector’s output is coupled with a Tektronix DPO4104 1.0 GHz digital oscilloscope 
which functions as an ADC and data collector, while a laptop computer running MATEAB 
performs the DSP. A ThorEabs EE502 filter is connected between the detector and the 
oscilloscope as an interpolation EPE with a 100 kHz cutoff frequency. All fiber used is 
SME28E single-mode optical fiber. 



Figure 5.1. DM-NYFR Flardware Setup with Significant Components Num¬ 
bered. The Components are Described in Table 5.1 
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Table 5.1. DM-NYFR Prototype Hardware Build Component Descriptions 


Number 

Part 

Name 

Function 

1 

S1FC11550 

ThorEabs CW laser 

1550 nm CW laser souree outputting at 1.75 mW. 

2 

DPO4104 

Tektronix seope 

Digital 1 GHz oseilloseope used as an ADC and data 
eolleetor when eonneeted to a laptop. 

3 

X5-MZI 

JDSU MZI 

Single-arm MZI biased at zero for pulse generation 
and positive quadrature for signal ineeption. The 
MZI has a 40 GHz bandwidth. One is eoupled with 
the EPG to generate sampling pulses, the other is 
eoupled with the RE signal input. Switehing voltage 

is 14 - 5.5U. 

4 

33220A 

Agilent funetion 
generator 

Programmable 20 MHz funetion generator. One is 

used to drive the EPG and another is used as an RE 

signal souree. 

5 

D-30IR 

Newport 

photodeteetor 

PIN photodiode with an internal 10 dB ENA that 
eonverts optieal signals into eleetrieal signals. 

6 

E3612A 

DC power supply 

Power supply used to power to the EPG eireuit. 

7 

EF502 

ThorEabs EPE 

100 kHz elliptieal low pass filter (5th order) used as 
an interpolation filter between the de tee tor and ADC. 

8 

EPG 

Eleetrieal pulse 
generator 

Generates eleetrieal pulses to drive MZIl. Design 

detailed in Seetion 5.3. 


5.2 MZI Setup and Testing 

The setup and testing of the MZIs used for modulating the CW laser is shown in 
this seetion. The JDSU X5 MZIs are mounted on silieon boards shown in Figure 5.2. Fiber 
attachments were soldered to the optical ports, and a coaxial cable attached to the RF input 
wire for the DC bias were attached to the MZI board. The MZI was mounted in a case to 
protect it from any physical damage. 
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Optical input 


Optical output 


Figure 5.2. JDSU X5 MZI Mounted on a PCB Showing the RF Input, DC 
Bias, Optical Input, and Optical Output 


The EPG circuit which outputs a 30.0-nanosecond wide electrical pulse at I 4 = 5.5V 
was connected to the MZI, and an output optical pulse was examined in order to show the 
MZI response. The oscilloscope capture of the MZI optical output is shown in Figure 

5.3. 
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Figure 5.3. Oscilloscope Capture of a Single 30 ns Wide Electrical Pulse 
Incident on a JDSU X5 MZI 
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The -3.0 dB pulse width was measured as 17.11 ns at 7.60 mV. The pulse-width variance 
was 303.1 ps, and the amplitude jitter was 126.0 fiY. The optical pulse width is smaller than 
the electrical pulse width due to the non-linear MZI input/output response (3.5). 

5.3 Pulse Generation Circuit 

An EPG circuit was built in the NPS photonics labratory using COTS components 
that can output a 30 ns wide electrical pulse train with an amplitude of 5.5V and a linear 
FM sweep as described in (3.8). The EPG is used to drive the first MZI and generate optical 
pulses sufficient for a prototype proof-of-concept build. 

5.3.1 Circuit Design and Simulation 

The EPG circuit was designed and simulated in NI Multisim, an electronic work¬ 
bench built around SPICE software simulation to accurately model circuit behavior. The 
schematic is shown in Figure 5.4 along with a graphical representation of how a FM sinusoid 
propagates through the circuit. The EPG consists of three main parts: a swept sinusoidal 
AC source, a comparator, and a one-shot multivibrator. 
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Figure 5.4. EPG Circuit Schematic Created in Nl Multisim Along with an 
Example Waveform 


Circuit Description 

The AC source is a chirp waveform output from a 33220A signal generator which 
is connected to a LT1016CN8 comparator referenced to ground. The AC signal is Vi and 
the reference level is V 2 . When Vi >¥2 the output is V^c, 5.5V. Otherwise, the signal is 
grounded to O.OV. This turns the chirp sine wave into a swept-square wave with distinct rising 
and falling edges. These edges are used as triggers for the pulse generator. The LT1016 
has unused reference pins that must be grounded in order to function correctly. The trigger 
signal from the comparator is sent to the SN74121N monostable one-shot multivibrator that 
outputs the electrical pulses. 

The SN74121N is a voltage oscillator that can flip between two states, Ycc and 
ground, when triggered by a clock edge. Only the ground state is stable, however, and the 
high state is unstable. A falling edge trigger pushes the multivibrator into the unstable state, 
setting the output high for a specific time determined by a programmable RC constant r. For 
the EPG, t=30 ns, and because the swept-square waveform functions as the trigger signal. 
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every falling edge outputs a 30 ns wide eleetrieal pulse with an amplitude of = 5.5V. 


NI Multisim Test 

The sehematie shown in Figure 5.4 was simulated in NI Multisim, and the results 
are shown in Figure 5.5. The SPICE models for the SN74121N and LT1016CN8 were 
downloaded and imported from an NI chipset database. The AC signal input and comparator 
output is shown with scope 1 and tests the square-wave conversion. The falling edge trigger 
from the LT1016 and a single 30 ns wide pulse is displayed by scope 2. The output has a 
perfect impulse shape due to the simulation parameters not including a realistic load. 




Scope 2 


Figure 5.5. Multisim Simulation Showing the Output of Scope One and Two 


5.3.2 Pulse Generation Testing 

The EPG circuit in Figure 5.4 was constructed on a breadboard for testing. The 
input was attached to a 1.0 MHz sine wave and the output connected to the test bench. An 
oscilloscope was used to probe various points of the circuit in order to test all the parts 
and confirm overall function. The oscilloscope outputs are shown in Figure 5.6. Shown in 
Figure 5.6(a) are probes connected to the input and output of the LT1016 converting the 
sinusoid into a clock signal. Additionally, the SN74121 functionality is demonstrated by 
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capturing the input clock signal and the output pulses, which are confirmed to be 30 ns long 
and are shown in Figure 5.6(b). 


(a) 



(b) 



Figure 5.6. Results of EPG Testing A) Clock Signal Generation with CHI 
Connected to the Input FM Sweep and CH2 Connected to the LT1016 Out¬ 
put and B) Pulse Generation with CHI Connected to SN74121 Input and 
CH2 Connected to the EPG Output 


A detailed examination of a single falling edge and timing pulse is shown in Figure 
5.7. These tests indicate the EPG circuit is functioning correctly. Lastly, the oscilloscope 
is attached to measure both the output electrical pulses and the output optical pulses. The 
MZI is coupled with a PIN photodiode in order to measure the optical pulses by converting 
them into electrical signals. This is shown in Figure 5.8 with both the electrical pulses and 
their corresponding optical pulses. The optical pulses slightly lag the electrical pulses, but 
otherwise, the MZI output response mimics the driving signal. 
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(CH2) 


The design, construetion, and testing of the hardware NYFR prototype and compo¬ 
nents was discussed in this chapter. The prototype used COTS components to implement a 
scaled down photonic NYFR that functions as a proof-of-concept build. Testing confirmed 
the individual components and their functionality. Test results from the functioning NYFR 
prototype with target RF signals is detailed in the next chapter. 
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CHAPTER 6: 

Prototype Test and Evaluation 


The hardware prototype testing, evaluation, and results of the DM-NYFR shown in 
Figure 5.1 is covered in this chapter. Multiple modulation schemes are tested with the first 
MZI stage to confirm the FM optical pulse train output. The primary results are from time- 
frequency spectrogram captures of target signals from the RF signal generator. Additionally, 
the signal processing methodology of de-noising, thresholding, and edge detection which 
made the incepted signals human readable is discussed in this chapter. 

6.1 FM Optical Pulse Testing 

Four FM presets are used as modulation schemes for the sampling pulse train. The 
FM presets are listed in Table 6.1 where tm is the modulation period, AF is the sweep 
bandwidth, and fc is the starting frequency for the sweep. All modulations are linear 
sawtooth chirps instead of triangular chirps, which were chosen due to hardware limitations 
with the signal generator. The varying bandwidths of the FM chirp signal presets allows 
a wide spectrum of target RF waveforms to be folded by the NYFR. Higher frequency 
signals are intercepted more efficiently with smaller sweep bandwidths, whereas larger 
sweep bandwidths perform better for lower frequency signals. The relationship between 
sweep bandwidth and compressed signal bandwidth is explored in detail in Section 6.3 
along with the effects of using a non-ideal modulation sweep bandwidth preset. 

Table 6.1. Pulse Train FM Presets 


Preset 

Parameters 

Preset 1 

im = 2 ms AF = 25 kHz fc = 100 kHz 

Preset 2 

= 2 ms AF = 15 kHz fc = 100 kHz 

Preset 3 

im = 2 ms AF = 10 kHz fc = 100 kHz 

Preset 4 

im = 2 ms AF = 5 kHz fc = 100 kHz 
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Oscilloscope captures are used to analyze the pulse-repetition interval (PRI) over 
a modulation sweep. This is shown in Figure 6.1, which highlights the changing PRI at 
the beginning of the modulation sweep and near the end of the modulation sweep. A PRI 
of 0.01 ms eorresponds to the starting frequeney of 100 kHz, while the PRI of 0.0081 ms 
eorresponds to a repetition frequency of 123.46 kHz. 


(a) (b) 






. ..id.. 

0.0081 ms 

< - > 



End of FM sweep: F 

= 123.46 kHz 


Figure 6.1. FM Optical Pulse Train PRI Measured at A) Beginning of the 
Sweep and B) Near the End of the Sweep 


Additionally, the captured oscilloscope output was analyzed in MATLAB in order 
to plot a time-frequency representation of the optical pulse train. The results are shown in 
Figure 6.2, whieh overlays the time-domain oseilloseope eapture with the frequeney-domain 
speetrogram to show the modulation bandwidth of 25 kHz. The test result confirmed the 
operation of the optical pulse generation cireuitry when the input signal is linearly chirped. 
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Figure 6.2. Time-Frequency Spectrogram of FM Optical Pulses Compared 
Against Amplitude-Time Oscilloscope Capture 



6.2 Signal Processing Methodology 

The signal processing methodology, specifically noise removal and thresholding, 
that is required due to the high noise floor of the optical signals is covered in this section. 
This is due to the low optical power on the input of MZI2 as well as characteristics associated 
with the MZI functionality. Each MZI has roughly a 3 dB loss between the input optical 
intensity and the output optical intensity. As the modulators are added in series, power 
loss increases, resulting in low incident power on the detectors. Another factor for the high 
noise floor when using MZIs is the imperfect extinction ratio that prevents total destructive 
interference. The ideal case of the MZI transfer function is illustrated in Figure 3.2; however. 
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in actual operation it is impossible to eompletely reduee the output intensity to zero and as a 
result some laser output is always present. Lastly, relative intensity noise (RIN) arises from 
power instabilities and vibrations within the CW laser eavity. With a low sampling signal 
duty eyele, the time-averaged noise power is larger than the signal power, and de-noising is 
a necessity. This is shown in Figure 6.3. 


■- Signal 


Optical noise 

Figure 6.3. Example NYFR Capture Illustrating Both the Signal and Optical 
Noise. 



The folded signal ineludes optical noise which needs to be filtered before analysis 
ean take plaee. Additionally, thresholding and edge deteetion are implemented to make the 
undersampled signals distinguishable. The image proeessing bloek diagram is shown in 
Figure 6.4 and details the DSP steps to refine the NYFR output. 



Figure 6.4. Block Diagram of NYFR Image Processing Steps 
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The raw ADC signal is first normalized and then passed into a noise removal bloek 
that uses squared averaging to eompute a eutoff value below whieh the signal is set to zero. 
The noise power Pnoise is defined as 


E'L, (s[<] - i)" 

^noise ~ 

n 


( 6 . 1 ) 


where 5[/] is the n sized signal array and i is the index eorresponding to a data point. The 
signal is de-noised by setting the signal array to 


r 



0 


if S[i] > P, 
if S[i] < P, 


( 6 . 2 ) 


whieh removes mueh of the optieal noise. After initial proeessing, the signal is thresholded 
in order to faeilitate edge deteetion. Otsu’s method of optimal thresholding is used to 
isolate the signal from the baekground noise by identifying the 0th and 1st order eumulative 
moments of an n-sized histogram. The maximum varianee of the elass separability between 
the two histogram moments is the optimal thresholding point [16]. This is shown graphieally 
in Figure 6.5, whieh is a histogram of a NYFR eapture with Otsu’s point labeled. Lastly, a 
simple edge deteetion method used binarization to eonvert the signal into a two-state array 
depending on the average power after thresholding. Anything less than the average power 
after thresholding was set to zero, whieh highlighted the folded signal slope. All signal 
proeessing steps were implemented with MATLAB. As a result of the signal proeessing 
steps, both the MZI optieal noise and RIN was effeetively removed from the ADC eaptures 
and the signals were able to be eorreetly unfolded. 


51 




Figure 6.5. Example NYFR Capture Flistogram Illustrating Otsu’s Method 
of Optimal Thresholding 


6.3 Single-Tone Test Results 

Selected DM-NYFR tests using a single frequency tone as the target RF signal 
are shown in this section. The sampling modulation presets used are listed in Table 6.1. 
Target RF signals are successfully undersampled and extracted using DSP, confirming the 
prototype’s functionality. 

6.3.1 Selected Tests 

Select tests with full DSP methodology are presented in this subsection in order to 
highlight the NYFR functionality. The first selected test was a 225 kHz RF signal with preset 
three as the modulation scheme. The results are shown in Figure 6.6, which details both the 
final signal and the image processing steps. All displays are time-frequency spectrograms 
computed in MATLAB showing a full AF frequency sweep and Nyquist zones A = 0, -1. 
The capture in Figure 6.6(a) is the raw received signal from the ADC block of the NYFR. 
The 225 kHz signal is masked by noise and is not discernible from the image. In Figure 
6 .6(b) we see the DM-NYFR signal after preliminary de-noising and normalization. At 
this stage the target signal is faintly visible from the background. The next capture. Figure 
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6 .6(c), is after the application of optimal thresholding, which removes almost all of the 
clutter and leaves the target signal clearly identifiable. Figure 6.6(d) is an illustration of the 
final output of the DSP block and is used in conjunction with the unfolding methodology 
described in Section 4.1 using (2.6) and (2.7). This is shown in Figure 6.7 as a red line fit to 
the slope of the signal, which was used to successfully extract the correct original frequency 
of the target RF waveform. 



(b) 




Figure 6.6. 225 kHz Test with Preset Three Showcasing the DSP Steps: A) 
is the Raw ADC Signal, B) is the Initial De-Nosing and Normalization, C) is 
the Optimal Thresholding, and D) is the Final Signal After Edge Detection. 
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Figure 6.7. 225 kHz Signal with a Fitted Slope for Unfolding 


The second selected test was a 525 kHz RF signal with preset four as the modulation 
scheme and is shown in Figure 6.8 as a series of time-frequency spectrograms. The results 
confirm the NYFR’s operation by replicating the 225 kHz test. Clear progression is shown 
through each DSP block culminating in a distinguishable signal that is unfolded into its 
original Nyquist zone. Interestingly, in Figure 6.7(a) we see the same weak tone signal 
in Nyquist zone N = -I that is present in Figure 6.6(a). This single frequency anomaly 
is not associated with any target signal and disappears when the noise is removed. This 
tone likely arises from a periodic oscillation in the optical noise stemming from the MZI. 
The modulators have a cosine-squared pattern transfer function as described in (3.3), which 
causes spurious frequencies. Another explanation is laser RIN vibrations resonating in the 
laser cavity. Although this does not appear to affect the DM-NYFR performance, these tones 
did not appear in the rSoft simulations, and this underscores the importance of hardware 
prototyping. 
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Figure 6.8. 525 kHz Test with Preset Four Showcasing the DSP Steps: A) 
is the Raw ADC Signal, B) is the Initial De-Nosing and Normalization, C) is 
the Optimal Thresholding, and D) is the Final Signal After Edge Detection. 


The modulation sweep bandwidth AF has a direet afFeet on the output modulation 
slope and ean adversely alter the results. A common occurrence when dealing with a large 
modulation slope is Nyquist zone overflow. The imparted modulation Sm is modeled as 


S„r = fc + MAF (6.3) 

where fc is the center frequency of modulation, AF is the sweep bandwidth, and M is a 
scale factor from the original signed Nyquist zone. As a result, grows as either the target 
RF frequency or modulation sweep increases. This can cause an overflow effect where a 
signal modulation crosses the bounds of a Nyquist zone, changing the slope value and sign. 
Depending on the capture size, we can detect and account for this overflow. This is shown in 
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Figure 6.9, where a 425 kHz signal is modulated with a 10 kHz AF. The overflow is evident 
at the 1.75 ms point, when the slope refleets back into the N = 0 Nyquist zone; however, 
due to the capture size of 4.0 ms, this was detected and corrected. It is recommended that 
the capture time F is at least twice the modulation time such that F > 2tm. 



0.5 1 1.5 2 2.5 3 3.5 

Time (ms) 

Figure 6.9. Nyquist Zone Overflow at 1.75 ms 


The hardware prototype testing and results for RF target signals was discussed in 
this chapter. The DM-NYFR functioned correctly, and applied signals were successfully 
unfolded into their correct original frequencies. This prototype, while not a fieldable 
design, serves as a proof-of-concept for photonic NYFR compressed sensing. The thesis 
conclusions and potential future work are discussed in the next chapter. 
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CHAPTER 7: 

Conclusion and Future Work 


The thesis results and future directions for further research are discussed in this 
chapter. Despite limitations on the physical prototype, the photonic NYFR architecture was 
successful, both in simulations and hardware, in sensing signals above of the traditionally 
required Nyquist rate. 

7.1 Conclusion 

In conclusion, this thesis research was successful in answering the question of 
whether integrated photonic components can be used to construct a novel compressed 
sensing receiver. Over the course of this thesis research, we implemented the NYFR 
architecture in the photonic domain, taking advantage of the increased bandwidths, sampling 
speeds, sensitivities, and EMF noise immunity that optical devices possess. Two potential 
design architectures, the DM and SM NYFR, were designed during the research. The SM- 
NYFR was primary envisioned as a technology demonstrator, while the DM-NYFR is the 
principal research contribution. The primary motivation for this thesis was to implement 
a novel undersampling technique with photonic components in order to take advantage of 
the large bandwidths of optical devices. This enables direct digitization and analysis of 
wideband signals from the RF spectrum. 

The design was successfully simulated in OPTSIM in order to confirm functionality. 
Computer simulations showcased the photonic NYFR’s ability to detect and recover signals 
an order of magnitude larger than the sampling frequency. The DM-NYFR uses two optical 
modulators. The first generates an FM optical impulse train, and the second incepts the 
signal as an envelope. The pulses are generated with an EPG outputting V„ amplitude 
30-ns wide pulses triggered every falling clock edge. The clock is frequency modulated 
with the specific NYER scheme. After the signals are modulated, they are coupled into a 
PIN photodiode sampled by an ADC at the Nyquist rate of 2cOs- Digital signal processing 
de-noised and thresholded the signals, which were then unfolded to recover their original 
Nyquist zone. Results with the hardware prototype show proof-of-concept and confirm 
the NYER architecture. Target signals tested included both single tone and pulsed RE 
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waveforms outside of the Nyquist eriterion, which were successfully unfolded. 


7.2 Future Work 

Suggested future work should focus on improving the hardware prototype beyond 
a proof-of-concept design. The primary limitations come from the COTS design of the 
pulse-generation mechanism. Described in Sections 3.4.2 and 3.4.3, the optical pulse 
width of MZIl is the biggest limiting factor that determines the NYFR upper interception 
frequency. Current EPGs can output high fidelity picosecond-wide pulses that can enable 
GHz-range target RF signals. Another pulse generation technique is to use mode-locked 
lasers to generate short optical pulses. Mode-locked lasers are commonly used to generate 
impulse trains; however, the time it takes for the laser to settle generally precludes frequency 
modulating the pulses. Novel PRF sweeping with mode-locked lasers should be explored 
further. 


Additional research directions include integrating the photonic NYFR with cur¬ 
rent photonic direction finding capabilities in order to develop a comprehensive wideband 
SIGINT/EFINT platform that would not only provide NYFR compressed sensing but also 
line-of-bearing (FOB) directional capability. Investigating a parallel channel configuration 
would allow multiple signal modulations per capture and would generate more informa¬ 
tion to unfold the signal. Using the NYFR as the front-end receiver, we could integrate 
with the wideband antenna system at F3 Technologies and implement autonomous signal 
classification and parameter extraction to test the photonic NYFRs capability in an inte¬ 
grated system. Fastly, implementing the photonic systems with superconducting quantum 
interference device (SQUID) modulators would improve the photonic NYFR sensitivity. 
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APPENDIX: MATLAB Code 


Nyquist Zone Code 
Richard Shmel, NFS 

computes the Nyquist zone from the slopes of the signal and modulation 


M = [0,-1,1,-2,2,-3,3,-4,4,-5,5,-6,6,-7,7,-8,8,-9,9,-IQ,1®,-11,11,-12,12]; 
Mf= [®,.5,1,1.5,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5,7,7.5,8,8.5,9,9.5,1®]; 

Mx = [®,-1,2,-3,4,-5,6,-7,8,-9,1®,-11,12,-13,14,-15,16,-17,18,-19,2®,-21]; 


rx notes 
rfxl = 4.594e-6; 
rfyl = l.®19e9; 
rfx2 = 5.9®6e-6; 
rfy2 = l.®61e9; 

S_clk = (rfy2-rfyl)/(rfx2-rfxl); %slope of the clock 
Signal 

sigxl = ®.4688e-6; 
sigyl = ®.1387e9; 
sigx2 = 1.688e-6; 
sigy2 = ®.368e9; 
center_freq = ®.15; 

S_sig = (sigy2-sigyl)/(sigx2-sigxl); %slope of the radar pulse (ie, rise/run) 

NYQUIST unfolding 

equations from section 2 

[X] = CS_sig-Mf"S_clk) .'^2 

[A,I] = min((S_sig-Mf*S_clk).^2); 

Nz = Mf(I); 

Nzf = Mf(I); 

Freq = Mf(I)+center_freq 
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Spectrum results 
Richard Shmel 

imports results and runs some basic de-nosing and thresholding 
described in section 6 


import file 
clear 

s1=’C:\Users\*\Desktop\PNYFR\v4\’; 
s2=’325k_5’; 
s3=’.csv’; 

s = strcatCsl,s2,s3); 

[X,Y] = importfile(s,14, 1QQ0013); 
fs = 16©e9; 

N = 150000; 
fs = 250e6; 

thresholding value from historgram 
thresh = hist_out; 
ax_l = 0.0; 
ax_h =0.1; 
load(’cmap.mat’); 

figure(l) 
hold on 
colormapChot) 

axis([-inf inf ax_l ax_h]) 

spectrogram(Y,hamming(N),9*N/10,N,fs,’yaxis’) 
savefigCstrcat (’ •'\Fold_test_l\ ’ , s2 , ’_raw’)) ; 
hold off 

norm = max(Y); 

Y = (Y./norm)* 5; 
new_max = max(Y); 

Y(Y<1.5)=0; 
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figure(2) 
hold on 
colormap(hot) 

axis([-inf inf ax_l ax_h]) 

spectrogram(Y.haitiitiingCN) , 9*N/10 ,N, fs, ’yaxis’) 
savefigCstrcat (’ ''\Fold_tGst_l\ ’ , s2 , ’_denoise ’)) ; 
hold off 

figure(3) 
hold on 
colormap(hot) 

axis([-inf inf ax_l ax_h]) 

spectrogram(Y,blackmanharris(N),9’'N/10,N,fs,’MinThreshold’.thresh,’yaxis’) 
savefig(strcat (’ •'\Fold_test_l\ ’ , s2 , ’_thresh’)) ; 
colorbar off; 
hold off 
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FMCW GEN 
Richard Shmel 

generates an FMCW sawtooth waveform to output to signal generator 
I&Q matrix 

can change bandwidth, period, and cycles 


sFMCW = phased.FMCWWaveformC’SweepBandwidth’,2Qe6,... 

’SampleRate’,l®0e6,’SweepDirection’,’up’,... 

’NumSweeps’,1,’SweepTime’,lQe-6); 
sig = step(sFMCW); 
windowlength = 16; 
noverlap = 64; 
nfft = 128; 

spectrogram(sig,windowlength,noverlap,nfft,sFMCW.SampleRate,’yaxis’) 
I = real(sig); 

Q = imag(sig); 

saveC’fmcwl.mat’,’I’,’Q’); 
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